clc
clear
options = odeset('RelTol',1e-2);
tic;
[T,Y] = ode15s(@testODEfunc, [0 10],[1 2 3], options);
toc;
% plot(T,Y(:,1))
tic;
[T,Y] = ode15s(@testODEfunc, [0:0.5 :10],[1 2 3], options);
toc;